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A bstract : 

An analysis of EPR line shapes by the method of Polnaszek, Brur.j, and 
Freed 1 h made for slowly tumbling vanadyl spin probes in viscous nematic 
liquid crystals. The use of typical vanadyl complexes as spin probes for 
nematic liquid crystals is shown to simplify the theoretical analysis und 
the subsequent interpretation. Rotational correlation times t and orienta- 
tional ordering parameters S where slow tumbling effects are expected to 
be observed In vanadyl EPR spectra are indicated in a plot. Spectral EPR 
line shapes were simulated for experimental spectra of vanadyl acetylaceto- 
natc (VOAA) in nematic butyl-p-(p-eLhoxyphenoxycarbonyl) phenyl carbonate 
(BEPC) and Phase V of KM Laboratories. Two factors are Important in achiev- 
ing good agreement with experimental spectra. The primary factor results 
from the a Low tumbling undergone by VOAA in these viscous solvencs. A 
secondary factor arises from the non-Brownian rotation of the small VOAA 
probe. Analysis of the inertial effects on the probe reorientation, which 
are induced by slowly flucuating torque components oi the local solvent 
structure, yield quantitative values for t and S, t . Similar effects were 
previously reported by Polnaszek and Freed for the small pcrdeuterated 2, 2,6,6- 
tetramethyl-4-piperldone N-oxide (PD-terapone) nliroxlde spin probe. The 
weakly ordered probe VOAA is in the slow tumbling region and displays these 
inertial effects throughout the nematic range of BEPC and Phase V. VOAA 
exhibits different reorientation behavior near the isotropic-nematic transition 
temperature (within r; l5°C) than that displayed far below this transition tem- 


perature. 
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Introduction 

The upplicut it'll of par iron Knot lc spin probe* in a variety of experiments 

has resulted In obtaining a wide range of structural and dynamical informa- 
1 2 

tlon. ’ Of particular interest to us is the use of complexes containing 
the vanadyl '.on VO" 4 as spin probes. * *’ In a recent paper, we had applied 
the stochastic Liouville method to calculate EPR line shapes in the slow 
tumbling region for vanadyl complexes in an isotropic medium.^ This anulysis 
allows one to extract quantitative Information from the F.PR line shapes for 
rotutlonal correlation times t in the range "10 - 10 ^ sec /rad. These 

slow tumbling EPR spectra occur when the vanadyl ion is rigidly attached to 
a macromolecule 1 or when a vanadyl complex is dissolved in a viscous 

i , H * 9 

solvent . 

We now want to extend our analysis of slowly tumbling vanadyl spin 

probes in an isotropic medium to that in an anisotropic medium. An example 

of current interest is the study by EPR of nem.it lc liquid crystals doped 

4,10 

with paramagnetic probes. ’ Appropriate analysis of the EPR spectra yields 

information of molecular orientational ordering and molecular reor ientational 
4 11-13 

dynamics. ' because liquid crystuls, especially the electro-optically 

14 

important room temperature nematogens, are highly viscous by nature, slow 

10 8 9 

tumbling EPR spectra have been obtained for both nitroxide* and vanadyl ’ 
spin probes. Since the slow tumbling region for vanadyl complexes occurs 
at a shorter T than for nitroxide spin probes, vanadyl complexes of 
comparable size to the nitroxide spin probe should more readily exhibit EPR 
spectra with slow tumbling features. 

A primary reason that vanadyl complexes are excellent probes for nematic 
liquid crystals is that they have nearly axially symmetric magnetic parameters 


4 


which la the same degree of symmetry possessed by the nematic liquid cryatala. 
The subsequent theoretical analysis is simplified and the results are easier 
to interpret than for typical nitroxide spin probes with less symmetry. 10 
Also, the characteristic eight hyperfine line EPR spectrum of vanadyl com- 
plexes provides a more severe testing of any theoretical analysis than a three 
hyperfine line nitroxide spectrum would provide. In addition, inhomogeneous 
broadening is usually neglible in the EPR spectra of vanadyl complexes. Pre- 
vious workers recognized these features of vanadyl spin probes but have en- 

countered difficulty in the quantitative interpretation of the EPR spectra 

8 9 

because of slow tumbling aspects. 

There are three main objectives for this paper. One is to establish 

12 1. 6 

when the breakdown of motionally narrowed formulae occurs. " * A second 

objective is to analyze the experimental vanadyl EPR line shapes by the 

stochastic Liouville method" as developed by Polnaszek, Bruno, and Freed 
13 

(PBF) for slow tumbling in an anisotropic liquid. Third is comparison 

of our vanadyl probe study to the detailed work of Polnaszek and Freed 10 who 

reported anamolous behavior for the perdeuterated 2 ,2,6,6-tet rameihyi-4-pi per- 

idone N-oxide (PD-tempone) nitroxide spin probe. 

Vanadyl acetylacetonate (VOAA) was used as the vanadyl spin probe 

because it has excellent EPR spectroscopic and physical properties and 

18 

because it has essentially axially symmetric magnetic parameters. An 
additional advantage arises from our previous detailed work with VOAA in a 
viscous isotropic medium. ^ The nematic mixture PhaBe V from EM Laboratories 
was studied because it is a nematic at low temperatures and has a wide 
nematic temperature range. A second low temperature nematic liquid crystal, 
butyl-p- (p-ethoxyphenoxycarbonyl) phenyl carbonate (BEPC), which Is a single 


4,15 
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component liquid crystal, was also studied for comparison. 

Experiment a 1_ 

Degassed samples of nematic liquid crystals doped with VOAA were 

prepared by standard vacuum line techniques. VOAA was purchased from 

19 

Alfa Products and was purified by rccrystal 1 lzat ion from acetone. The 
iicm.it lc liquid crystal HKPC was obtained from Kastman Cnemicals and was 
recrystallized from methanol and then from benzene.* 1 ' The nematic mixture 
Phase V was purchased from Ut Laboratories and used without further purlfi- 
cat ion. 

EPR spectra were taken on a Varian V-4500 spectrometer with 100 kHz 
modulation. Varian K-257 variable temperature unit controlled the tempera- 
ture to ‘1 K over the active region of the cavity. The magnetic field was 
measured with a Brucker BNM-12 tracking magnetometer, which was used for 
spectral calibration. The microwave frequency was obtained with a Hewlett 

Packard 5248 L electronic counter and 5257A transfer oscillator. 

Theory 

A. Slow Tumbling 
13 

PBF have developed a general method for calculating slow tumbling 

EPR line shapes in an anisotropic liquid. This approach ij based upon the 

stochastic Liouville method*^’" and is a generalization of an analysis by Nordio and 

co-vx>rkirs for the motional ly narrowed EPR line wldtlisof a Brownian particle 

12 

in a liquid crystal. Here, we apply the general PBF approach to vanadyl 
spin probes in viscous nematic liquid crystals. 

Previous detailed studies in an isotropic solvent have shown that the 
spin Hamiltonian for VOAA has very nearly axially symmetric magnetic inter- 
actions.^’^ The mot tonally narrowed line widths (x < 8 x 10 " sec/rad) 


b 


calculated with axially symmetrii parameters differ by less than 0.5Z trom 

18 

those calculated with completely asymmetric trnsir values. Our slow tum- 
bling study of VOAA In toluene showed that the use of axially symmetric 

parameters gave very good agreement with experimental line shapes for the 

-11 -0 . 7 

slow tumbling range of t * 8 X 10 - 1 x 10 sec/rad. Only for EPK 

*8 

line shapes (t > 1 x 10 sec/rad) very close to the rigid limit spectrum 
would the use of completely asvmmetric parameters improve the fit to ex- 

. i 21 

perimental spectra. 

The orientational portion of the pseudo-potential u(0) for nematic 
liquid crystals may be expressed by an expansion 


u(o) - 2 _, \ C08 " 0 (1) 

n even 

22 

where 0 is the angle between tlie molecular long axis and the director. 

Because VOAA hat the same degree of symmetry as the nematic liquid crystal- 
line solvent, the form of Eq. (1) should also reflect the orientational restoring 
potential for VOAA. In our experiments, the small VOAA is moderately ordered. 

Previous studies have shown that El’R spectral line shapes are not sensitive 

13 23 

to variations. Thus a simple Maier-Saupe potential 


u' (0) • Y 2 cos'O 


( 2 ) 


is suitable to describe VOAA. (The prime refers to the spin probe values.) 
In Eq. (2), 6 is the angle between the VO bond (molecular z axis) and the 
director. Note that VOAA, an oblate symmetric top molecule, has 0 ■ 

so thut Y is positive. Eq. (2) leads to a Boltzmann equilibrium P^(0) 


2 


7 


given by 


12 


j 


| J 

P > (o ) • exp (X com 0)/ do Bin© exp (X cob 6) 


where 


X - -y;/kT . X < 0 for VOAa. 


(3) 


(4) 


23 

An orientntion.il ordering parameter is then defined by 


S * m < r 2 <ll) > 


(9 


-r 

«' o 


dOsinO P„ (6) p (0) 
l o 


where 


( 6 ) 


P, (9) - Y <3 cob 2 « * 1) 


(7) 


71 i 

For complete alignment where P y i©) - / (0 - —) , • - — . 

For a Brovnlen particle under a restoring potential Kq. (2), matrix 
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elements of the symmetric rotational diffusion operator T are given by 

N(L,L' ) L I - 1 L' x. 2 2 2 

— 2“< D fl. r V >- lR l L(L * l) + (R 11 ‘ R l )n + I5 V 1 6 LL’ 6 mm* 

8n 1 ' 


- 2R X »(1 - JjX-l)* N(L.L')| „ l .1 o I’j 


ram 


(-!)■ Ha.O)(ooo') (-mO») 


5 , 

mm 


( 8 ) 


Eq. (8) is appropriate for anisotropic viscosity where R^and j are rota- 
tional diffusion constants perpendicular and parallel, respectively, to the 

1/2 L 

director. Normalization factor N(L L') is [ (2L + 1)(2L' + 1)] D are 

om 


a 


Wlgner rotation matrices, and ^ j are 3-J symbols." 4 Note that effect* 
of anisotropic viscosity are proportional to (R^ - K^m*. 

The application of the general methods In PBF to the VOM spin probe 
Indicate It* simplified nature. In addition to the simple one-term potential 
Eq. (2), VO M has the principal axes of Its magnetic parameters (g and A 

tensors) coinciding with the molecular rotational coordinate frame. One 
need not Incorporate tilt angles and additional transformat i- as of coordinate 
syms terns into the orientation-dependent portion of the spin Hamiltonian 
(other than usual one that is doiu even for Isotropic solvents). 

Appendix A explicitly presents a convenient set of equations whose 
solution yields the unsaturated EPR line shapes for vanadyl spin probes 
undergoing Brownian rotation (Eq.8 ) in nematic liquid crystal. Nonsecular 
contributions, in addition to those given for the laotropic solvent, are 
also Included. In these equations, the director Is parallel to the external 

25 

magnetic field. Tills alignment is the usual case in nematic liquid crystals. 

A second method of solution allows tor inodif ication** that arise from 

non-Brownian diffusion. This method Involves the numerical diagonalization 

of Eq. 8 (a real symmetric matrix) to obtain eigenvalues E n and eigenfunctions 

os 

n 

tt • 

om 


-y 

om /_> 


a (l,)L l) L 
om otn 


(9a) 


E n - < u n r u n N 
ora \ om om / 


(9b) 


(n)L 


where the coefficients a form 

om 


column vectors of the post-matrix 


in the diagonalization transformation. In general, eight diagonal izat ions 
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would be performed for vanadyl (corresponding to the eight m values). Devia- 
tions from Brownian rotation may be expected for the small VOAA probe dissolved 

In the much large! liquid crystal solvent moleruleH. One worthwhile modifi- 

26 

cation is an analysis of inertial effects" to lowest order, which are in- 
duced on VOAA by sh>wly flucuutlng torque components from the local solvent 
27 28 10 

structure." ’" Polnaszek and Freed applied this analysis to nitroxlde 

29 

spin probes embedded in a liquid crystal solvent. Polnaszek gives details 

for obtaining equations for the s?ow tumbling region. Specific equations 

are given in Appendix B. An essential point of this approach Is that factors 

c' (> 1) are introduced to Indie ite deviations from Brownian rotation. 

L, m 

Values i | m ■ 1 yield Brownian rotation diffusion. In a motional ly narrowed 

analysis, t' would correspond to correction factors in the spectral density 

£ ,m 

functions. For K. - K,, ■ and low ordering, c' _ (r' ) and e' , n (t' ) 

JL 11 bt 2,0 sec 2,mf0 ' psec 

reflect the relative effective torque components (i.e., parallel and perpendi- 
cular, respectively, to the director) which affect reorientation of the probe. 

B. Motionally Narrowed Region 
13 

PBF describe a convenient method for obtaining line widths from the real 
part of the spectral density function K(m,w) given by 


He K(m,u) 


V ( ' 1)Be - - 
'Vo e 2 + „ 2<U » 

nrO n 


om 


u ><u 

n 


D 


o,-m 


u > 


( 10 ) 


where and u^ are the eigenvalues and eigenfunctions of Eq. (9). Those non- 
secular terms which arise from the large intramolecular magnetic interactions 


are treated firstly by a perturbation analysis. 


7,13 


Linewidth contributions 
13 

from these terms are then obtained by the method in PBF. Such an analysis 
yields a cubic dependency (6n/) in the line width expression/ 4 ’ ‘ ' * H The 
first derivative line width (in gauss), which is slmilur to the one given 


for isotropic solvents, may be written 


A H (m) 


a + 8® + ym + 6m" 


( 11 ) 
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where 


- o" ♦ 


a ■ a 


l Ke 


r\y\ 


V' [K(0,0) ♦ l K(0,ui ) ♦ - K(0,0)J 

4 o 14 u 


+ 1(1 + 1)C 2 [K(1,» ) 2K (2,u ) + ~ K(0,w ) 
a o 3 o 

+ /I ffln (K(0i0) . 1 KH, U> )]| 


./31 r 


3 KU[K(0,0) + | F(l,u» o ) + ~ K<0.0) ] 


n IyJ 

+ A l t 1 . _ + -Uj JL [k(1,ui ) - ~ K(0,0)J 

w a J 

+ -- Id + 1) — l K(0,0) - K(l,w) J 

, 5 J 2 I) 3 , Da v/>1 . 

7 ) iii a w a 


Y “ 


2 Ke 


/3 |y | 
1 «• 


3-lK(0,0) - j K(l,u. a ) 


i K(0,w ) 

o o 


+ I K(1 ' w o ) “ 4 K(2 * w o ) 




- ~ IK(0,0) - | K(l.u» a )] 


2 Re 


IyJ • 


2 

4 [K(0,0) - \ K(1 »w ) ] 

3 a) 4 a 

o 


- fx 7[<(W)-§<o,.,)i 

o 


(12a) 


(12b) 


(12c) 


(1 2d) 


F.D.a.ci" are defined as in Kef. 7; u ■ 4 a; and u is defintd here as the 

II O 

klystron frequency. 

S in Kq. (5) is obtained readily by* h 
2 


s 

z 



(13) 


where <n> is the observed separation between hyperfine lines, A, ^ is the z 
principal A tensor component, and a' is the isotropic splitting corrected 
for nonsecular shifts. 

These shifts (in gauss) are given by the formula 


A B(m,S z ) 


A B(m) - 


S 

7 


U) 

X 


nU 





id + i) 


+ /| aDI(I + 1) - Dm - rfm J 

- /\ aDm 2 \ 


( 14 ) 


where A B (m) is the shift in an isotropic solvent. 
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Results and Dis cussion 

A. Magnetic Parameters 

The extraction of Infer t Ion about molecular ordering and rotational 
dlffualon from KPK line shapes requires that the components of the g and A 
tensors In the spin Hamiltonian he known. These magnetic parameters were 
determined from the combination of an isotropic rigid spectrum and an iso- 
tropic mot tonally-narrowed spectrum. The rigid isotropic spectrum was 
obtained by cooling the liquid crystal sample through the nematic range with 
the magnetic field turned off and then by recording the spectrum at 77K. 

Aj j and g ( j components are readily determined from this spectrum. The iso- 
tropic motional ly narrowed spectrum was taken at 110®C, which is well above 
the nematic-isotropic nematic transition temperature. With the a and g ( 
values determined from this spectrum, A L and g a components are gotten from 
the equations A^_ - — (3a - A^) and g x • j (3g (> - gjj). The magnetic para- 
meters for VOA A were the same in both BEPC and Phase V within experimental 

error. The values are Ajj ■ -185 i 2 C, ■ - 68 1 1C, gj ^ • 1.943 t .001, 

18 

and g^ - 1.982 ♦ .001 and are essentially those obtained in toluene. 

B. Motlonally Narrowed Region 

The approximate region where a raotlonally-narrowed analysis is applicable 
for vanadyl spin probes Is Indicated In Figure 1. Figure 1 Is a theoretical 
plot of t versus S z where the motlonally narrowed region Is the area below 
the appropriate line. The solid line represents the breakdown of the motion- 
ally narrowed F.q. (13) for the determination of the orientational ordering 
parameter S^. Only negative values for S were considered since typical 
vanadyl spin probes have 8 ■ fl/2 (VO bond perpendicular to the director). 

The dashed line Indicates the breakdown of motlonally narrowed line width 


n 


formula Eq. (11) for obtaining rotational dynamics. Determination of these 
breakdowns was made by comparison of motionally narrowed results to the 
general slow tumbling results (Appendix A), which are also valid in the 
motionally narrowed region. Calculation were for isotropic Brownian rota- 
tional diffusion (R, . * R. = -r- ). The criteria for the lines dr . v.i in 
11 i. bT 

Figure 1 are IX deviations in the line separations and 3% deviations in the 
line widths from the motionally narrowed calculations. For a low degree 
of ordering, these deviations begin to occur in the high-field outside 
hyperfine lines. For high degree of ordering, deviations occur initially 
in the central lines. VOAA magnetic parameters were used in these calcu- 
lations. Slight differences in Figure 1 would result from using 
magnetic parameters of other vanadyl spin probes. Also, rotational diffusion 
other than the isotropic Brownian model would show slight discrepencies from 
that given in Figure 1. The use of Figure 1 presumes a recognition of the 
sensitivity of EPR spectral line shapes to changes in t and i(S z ). When the 
vanadyl spin probe is highly ordered, the EPR spectrum will naturally become 
less sensitive to rotational diffusion and ordering potential changes. Such 
changes, in practice, may be camouflaged by an orientation-independent line 
width a" which may be regarded as a resolution parameter. 

C. Slow Tumbling EPR Line Shapes 

Experimental EPR spectra for VOAA were recorded in the nematic range 
from 82°C to 37°C for BEPC and from 72°C to -29°C for Phase V. Figures 2-4 
show typical experimental spectra that are both near the isotropic-nematic 
transition temperature T^ and reasonably near to the nematic-solid transition 
temperature. Simulated spectra calculated from equations in either Appendix 
A or Appendix B arc presented as broken lines. All experimental spectra 
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formula Eq. (11) for obtaining rotational dynamics. Determination of these 
breakdowns was made by comparison of motional ly narrowed results to tiie 
general slow tumbling results (Appendix A), which are also valid in the 
nationally narrowed region. Calculation were for Isotropic Brownian rota- 
tional diffusion (Rjj • ). The criteria for the lines drawn in 

Figure 1 are 1% deviations in the line separations and 3Z deviations In the 
line widths from the motlonully narrowed calculations. For a low degree 
of ordering, these deviations begin to occur In the high-flcld outside 
hyperflne lines. For high degree of ordering, deviations occur Initially 
In the central lines. VOAA magnetic parameters were used in these calcu- 
lations. Slight differences in Figure 1 would result from using 
magnetic parameters of other vanadyl spin probes. Also, rotational diffusion 
other than the isotropic Brownian model would Hhow slight dlscrc pencics from 
that given in Fi l. The use of Figure 1 presumes a recognition of the 

sensitivity oi EPK spectral line shapes to changes in T and A(S ). When the 

z 

vanadyl spin probe Is highly ordered, the EPR spectrum will naturally become 
less sensitive to rotational diffusion and ordering potential changes. Such 
changes, In practice, may be camouflaged by an orientation-independent line 
width a" (T,, ^ ), which may be regarded as a resolution parameter. 

C. Slow Tumbling EPR Line Shapes 

Experimental EPR spectra for VOAA were recorded in the nematic range 
from 82®C to 37°C for BEPC and from 72®C to -29®C for Phase V. Figures 2-4 
show typical experimental spectra that are both near the isotropic-nematic 
transition temperature T and reasonably near to the nematic-solid transition 
temperature. Simulated spectra calculated from equations in either Appendix 
A or Appendix B are presented as broken lines. All experimental spectra 
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taken In the nematic range clearly exhibited the eight hyperfine line spec- 
trum indicative of motional averaging. These line shapes suggest that there 
should be enough rotational motion about the VO bond to average out the 
small magnetic anisotropy in the x and y principal axes. Subsequent 
analysis bears this out. Thus the u.'c of axially symmetric magnetic parameters 
and of F.q. (2) for the restoring potential is feit to be well justified. All 
spectra In the nematic range had highly asymmetric lines which ure very sug- 
gestive of slow tumbling effects. In fact, spectra In the isotropic phase Just 
above the isotropic-nematic transition temperature T^ were obviously in the 
slow tumbling region. Although there was considerably supercooling below 
the nematic-solid freezing point, no discontinuities were observed in the 
EPK spectra or in their subsequent analysis. Both the observed isotropic- 
nematic and observed nematic-sol Id transition temperatures were clearly 
discerned by dramatic changes is the KPR spectra. 

The computer simulations for the slow tumbling EPR line shapes are 

13 

calculated by the method developed by PBF (See Appendix A for the explicit 
equations that were used) or by a modified vers Ion ^ U (See Appendix B). 

Fiquire 2 allows the experimental spectrum for VOAA in BEPC at 79°C, 

which is Just below to T. . Simulated spectra, that are calculated for Brownian 

rotational diffusion (Appendix A), are also shown. Spectra for Isotropic 

Brownian rotation are simulated by varying essentially only two parameters, 

T and A. The general effects of varying T and A on spectral line shapes 
13 

are given in PBF . Note that for vanadyl spin probes, which have A < 0 and 
|A^J < |Ajj |, the effects of slow tumbling is to shift intensity of the 
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hyperflne lines c loner together. Thun a not tonally narrowed analysis, which 
would ascribe such shifts to an Increased ordering (X more negative) rather 
than to a longer t, overestimates the degree of ordering for the slow tumbling 
region. Figure 2 Ik u typical best fit of an isotropic Brownian rotational 
diffusion model to an experimental VOAA EPR spectrum. Although significant 
dlscrapuncles arc Indicated, the out invited best fit S and x values 
are reasonable and represent a significant Improvement over i motlonally 
narrowed analysis. Improved fits to the experimental spectra were made by 
Introducing *n anisotropic viscosity (R^ i Rj). Spectral simulation for 
thic model were fit by varying essentially three parameters (t^, x^, and A). 
Figure 2 shows the improved agreement with x^ » Hx^. However, the ratio 

Increases rapidly and x^ becomes shorter with decreasing temperatures. 
Similar implausible results for nltroxidc probes have been reported by 
Polnaszck and Freed,' 0 who then invoked a slowly flucuiting torque analysis. 

Such an analysis is expected to be appropriate here because of the small 
size of VOAA relative to the ordered liquid crystalline solvating molecules. 
Figures 3 and 4 show spectral fits which are based upon this analysis (Appen- 
dix B). For purposes of comparison to nltroxlde spin probe studies, only 

two correction factors, c’ and c' , to spectral density functions will 

sec pset 

be considered. Note l’ ■ l yields Brownian diffusion results while c' > 1 
indicates deviations from Brownian diffusion. Also we let R^(0) ■ Rj (0) ■ -~ 
The intention of this analysis Is to establish a set. of c' values for all 
temperatures In the nematic range. With c' values fixed, spectral fits arc 
made by varying two parameters, x and A. The values ■ 1.5 t 0.5 and 

l pse . “ 10 t 2 wuu ld appear to be a consistent set of parameters except for 
a *-15°C temperature range below the T^, • (See Figures 3 and 4.) For 
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this temperature rang*, 1 p Sl , ( ■ 4 t 2 gave better fits. Although aome 
systematic dlicrrponcleH between calculated and exper Inent il spectra h 1 1 1 1 
exist, the spectral fits were good. Tills analysis gave 1 values that were 
precise to 101 except for low leuiperutures In Phase V where the precision 

was 202. S values computed numerically from V (See Eq. 3) were precise 

Z 

to *52. 

The 1 and S values determined from these best fitting spectra, which 

Incorporated the i' corrections, are used In Figures 5 and 6. Figure 5 is 

a plot of S versus temperature (T) where S is calculated numerically from 
z z 

Eq. 3. A plot of log T versus ,j. Is shown In Figure 6. With the neglect 
of a 13°C temperature range below T^, Figure 6 for Phase V shows good Ar- 
rhenius behavior with an activation energy of 9.0 ♦ 0.4 kcal/mole. This value 

10 29 30 

is comparable to that found for nltroxldc spin probes. * * The pre- 

exponential factor is (3.4il.l)xl0 * f> sec/rad. BEPC had too short of a 
nematic range to determine accurate Arrhenius parameters. 

Attempts to account for the deviation from Brownian rotation by in- 
voking a static distribution of director orientations were not successful. 

Here, i line shape g(Aio,| ) was calculated by the aummatlor (i.e., a trape- 
zoidal rule) of spectral amplitudes f(Au),<J>) weighed by the distribution 
2 2 

factor exp(-sin <$/t^)sl*i where is the angle between the director and 

10 32 

the external magnetic field. * $ is the root mean square deviation (spread) 

in $ . In particular, the experimental spectrum for VOAA in BEPC at 79 C 

was used. For Brownian rotation with T ■ 4.1 x 10 *°sec/rad, ■ -0.13, 
and 4> ^ ■ 0.17 rad (10 ); the calculated line shupe g(Au),4> 0 ) showed a rhange 
in asymmetry of the hyperf ine lines opposite to that which Is observed. 
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Director f lucuat ion* are predicted to have ncgllble effects on Llie 
LPK spectra of the rc.ulively small VOAA probe. I( * Also, a simple Jump model 
of rotational reorientation was inadequate to explain the discrepencies from 
brownian rotational diffusion. 

Conclusion 

This work Is an Initial Inquiry of slowly tumbling vanadyl spin probes 

In an anisotropic medium. Here the stochastic Llouville method as developed 
1 3 

by PBF has been appropriately applied to the analysis of EPR spectra from 
VOAA In viscous nematic liquid crystals. The use of vanadyl complexes as 
spin probes for nematic liquid crystals greatly simplifies theoretical cal- 
culations and interpretation. Oulde lines for the breakdown of motlonally 
narrowed formulas were established (see Figure 1). 

Spectral simulations for experimental spectra of VOAA In nematic BEPC 
and Phase V showed that VOAA was in the slow tumbling region throughout the 
nematic range of both these viscous solvents. Also, deviations from Brownian 
rotation were noted. An analysis, in which inertial effects on the probe 
reorier nation are induced by slowly ilucuating torque components of the local 
solvent structure, results in improved spectral fits. Sich an analysis, 

which introduces correction factors c' und t' , yields quantitative values 

sec psec 

for t and S^. t values for the weakly ordered VOAA showed relatively little 

variation with temperature near the isotropic-nematic transition temperature 

T„. (See Figure 6.) Here, c' ■ 1.5 t 0.5 and c' ■ 4 ± 2 gave good 

spectral fits. For temperatures more than 15°C below T^, x showed good 

Arrhenius behavior with an activation energy of J.O t 0-4 kcal/mole in Phase 

V. Here e' - 1.3 ± 0.5 and e' ■ 10 1 2 gave good spectral fits. These 
sec psec 


1H 


results suggest that nomat often* undergo orientational flucuatlons over a 
reasonably wide temperature range near that are different from the fluctua- 
tions that may occur tnr below T^. Our studies Indicate that vanadyl com- 
plexes may be good probes tor the study of phenomena near 1'. Specifically, 
VOAA was obviously in the slow tumbling region, where EPR spectra are most 
sensitive to motional dynamics, at of BEPC and Phase V. 

In their detailed study with the small nltroxide spin probe PD-tompone, 
Polnasxek and Freed* Introduced a slowly flucuatlng torque analysis to 

account for "anomoious" EPR line shapes. In particular, values of i' ■ 1-2 

sec 

and ■ 15“20 were used In the Incipient slow tumbling region (T < -6"C) 

In *hase V. Results for VOAA in this region show similar effects but with 

smaller t' values. This difference may be attributed to a larger effective 

rotational radius tor VOAA. However, our results indicate slowly flucuatlng 

torque effects throughout the nematic range but much reduced near T . Other 

K 

comparisons show that VOAA exhibits behavior similar to that of PD-tompone . 

V JAA Is sllghlv more ordered and It reorients more slowly with roughly the 

. . 10 

same activation energy. 

In summary, our results Indicate that both the slow tumbling aspects and 
the effects of non-Brownlan rotation should be resolved in order to extract 
quantitative Information about molecular ordering and rotational mobility. 

This present study should not be considered a critique of the slowly flucuatlng 
torque analysis. The work here does show that, effects which wore observed 
tor PD-tempone, are also seen for the commonly used probe VOAA. Undoubtably, 
the slowly flucuatlng torque analysis does provide a proper framework for 
explaining such effects. Further work is needed to characterize the correction 
f ailor c' (e.g. different size and shape probes). Refined theory and develop- 


ment of ocher techniques of solution should yield a clearer and more precise 

explanation for the motional dynamics o f the spin probe. Hucause of the 

10 29 

present necessary simplification In theory, ' inclusion of other effects 
(e.g., director distribution and flucuatlon, inlsotroplc viscosity) into 
the slowly flucuatlng torque analysis were not Justified. 
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App endix A 

This appendix indicates explicitly the equations that were used for 

obtaining vanadyl EPR line shapes for Brownian rotational diffusion with an 

misotroplc viscosity (Eq. H) under a simple one term Maler-Saupe restoring 

potential (Eq. 2). These equations may be obtained most easily by modlfica- 

13 

tion of the Isotropic liquid equations, which are explicitly given by 
Eq. A.l - A. 4, B. 3 of Ref. 7. The following term is added to the left-hand 
side of Eq. A.l - A.l in Ref. 7: 


-U(R U - ».)■>• ♦ T3 K. ‘ 2 1<„0> 


♦ U R, *U - ><-»■£ N<L.L')( o „ Jj') / l „ L ') <£„ 

L* V ' 


(J) 


+ 1 if R 1 X2 <- 1)P 


1 0 0 o') (-iol') C «> < A ' l > 

L' ' 


The right hand side of Eq. A.l of Ref. 7 becomes 


u ■ [~p— I)' - (o, G,o)exp(*r A cos’ ')sinl'd'J (A. 2) 

l- I / 0,0 z 

o J „ 


where I i3 

o 


'I 


exp (-^ A cos^ 6) sin Odd 


(A. 3) 


The absorption Z" In Eq. (A. 4) of Ref. 7 now becomes 


Z" • la Y ^ u,C L (2m +21 + 2, 2ra + 21 + 1) 
Zj l- °»° 


L m*-I 


(A. 4) 


2 ! 


Contributions from nonscculir forms In addition to those given by expression 
B.3 in Hot'. 7 have boon Included. Here, nonsecular corrections to the "for- 
bidden" transition terms have been made. Nousccuiar contributions to the 
coupling cjr n (i) • c!^ "* (i) for allowed transitions, which are made significant 

U|U V| Wj 

by Eq. (A.l), have also been included. These nonsccular contributions are 

7,13 

obtained by straightforward application of previous perturbation analysis.' 

The results for the coupling to the c|‘ (a,b;c,d) transition terms are given 

u j n 

by the expression 
. n+1 


(-O' 


1 ) 3F 2 /6(m+n)FD . 1>‘C 141(1+1) - 2m 2 - 2nm - n 2 ] 

w |_ | 20 " 5 30 


♦ 21(1+1) - 2m“ - 2nm - n“ J ^ c|j ^(a.bjc.d) 


+ 16n.‘ + 16nm 


+ £ N(L.L') - *y (m + f )F ° + ^[*-101(1+1) 

+ 8nJ + ;‘|L-2HI+1) + 2 m J + 2 nm + n^||j 2 2 Jj j^o,n (a ' b;C 

+ N(I.,L') J- 

L * l 


,d) 


+ + 2 )KD + 121(1+1) - 36m 


- 36nm - 18n' 


V 


i} ( s s s*)(-i 4 o l :) s>*.. 


d) 


+ ^ N(L,L')f(l,m+n) 


L' 


[- 


5Kb ^6(10m + lOn + 5)D 

28 84 


_baF 
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(2m + 2n + l)aD~| L2L'l L 2 L' UL' , . . , , 

1 — r — J | o o o -,.i »*i | Vi ( * tl ' klc ' d,,) 

S' .fi^OFl) V5(2n + 2n + 1)D 2 ] 

L 70 33 -1 


+ ^ N(L,L')f(I,i*+n) 

L' 

'l 4 L’ 
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3* 30 FD 3»5(2m -l)D"l L 4 L'l L 4 L'UL' , . . , ,v 

- TT J|o 0 0 -n -* -'l C « - , ^ ;c-l ,d) 


1 irfl j 0,n+l 


^ M/1 |tw/f , v f 5FD /6(10m +10n - 5)D 2 . •'baF 

(2m +2n - l)aDl / L 2 L' 1 I L 2 L' xL' , 
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la 
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(A. 5) 


The above expression (A. 5) Is added to the left hand side of Eq.(A.3) in 

Ref. 7. Eq.(A.2) In Ref. 7 may be regarded as a special case of Eq.(A.3) 

in Ref. 7. Expression (A. 5) was also used for the additional nonsccular 

coupling to the allowed transition terms. In the derivation of expression (A. 5), 

R << u>. Also, within second order perturbation theory, the nonsecular coupling 
_1 • _ 

couplings e' , (k) between C (i) and C_ ,(k) were approximated by 
w»n u, ii u f n 

the equations 


-L’ 1, L' L' . 

’ 0,n “ 2 ( L 0.n ' 0,-n ) 


(A. 6a) 


“ L ' 1/ L' 

L 0,nll “ 2^ t 0,nil 


L’ 1 

C 0,~(nilK 


Note that for n-0. C 0#0 -/2C 0#0 and - /2 e Q j • 


(A. 6b) 


To humw.ir Ize, the equations which wire used are obtained from Eqs. (A.l-A.3) 

in Ref. 7 that have been modified by (B.3) of Ref. 7, (A.l-A.3), and (A. 5). 

Eq. (A. 4) yields the absorption line shape. Kqs. (A. 2) and (A. 3) are 

12 

solved numerically with recursion relations. The first derivative EPR 

. 3 

line shapes are then most efficiently obtained by a diagonal lzation method. 
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Appendix Ji 

Details for explicitly obtaining slowly f Imitating torque effect* In 

29 

the EPR line shapes In given by Polnnnzek. I'lie general expression way 
be written 


(w - u - lq t" 1 ) C" (j) + <u" ! [H. (ft) , c”' , | u" ' , > 

J ^ q o,m,p ,q J , ^ o,.n'‘ 1 o.ra'.p.q' o.m' 


n ,m 


1/2 


( rTTc n , <j) + ^c n , (J)l 

\ R I / P o,m,p+l ,q J m,p-l,q 


♦ 



(*P 


/q + 1 C 


n 

o 1 ,q + l 


(J ) + ' P ♦ 1 *'q C" . . ,(1)1 

J 1 1 o,m,p+l,q-l VJ/ * 


■ 6 

no 


where t 

q 


6 « 6 
mo po qo 


//' - 1 
6R 


B.l 


B. 2 


and represents a correlation time for the decay of the nertlal effects 
Induced by the solvent. To Include effects Indued by torque components 
parallel and perpendicular to a director of a nematic liquid crystal solvent, 
c' may be labeled l ^... a nd l p. . to express these parallel and perpendicular 
components, respectively. Note that for i ' 1, Brownian rotational diffusion 

Is obtained, t' > 1 Indicate deviations from Brownian rotation. 

Eq. (B.l) has only p ■ 0, q • 0; p ■ l t q ■ 0; and p ■ 0, q » 1 combina- 
tions for the lowest order inertial effects. Note that the eigenfunctions 

u" and eigenvalues E° are used (Eq. 9). 
o,ra o,m 

The procedure to obtain the equations used for the slowly flucuatlng 

torque analysis may be summarized in three steps. First, diagonalize the 

rotational diffusion operator (Eq. 8) to obtain u~ and E*! .(Sec Eq. 9.) 

U,m U, m 


The real symmetric rotational diffusion matrix (Kq. 8) t4ay bo blocked according 
to different m values and each block may be numorlcnlly dlngonallxed sepa- 
rately. Second, pcrfoim n similarity traits format ion on the matrix formed 
from the coefficients of the transition terms cf' developed In Appund<x A. 
Retain only those terms that have F or D in them. The columns of the 
similarity transformation post-matrix are the vectors whose components are 
the coefficients a,V * in Eq. 9a. This step calculates the summation (second 

U | HI 

term) on the left hand side of Eq. D.l. Third, form the entire matrix by 

using Eq. H. 1 . This matrix may he viewed as partitioned Into thre sections 

according to the three combinations of p,q that were used. 

The absorption line »hape Z" Is then given by the following 
* 

Z" - Ira V cj’ 0 0 0 (2nrf2I+2,2m+2I+l) . B.3 

m »- 1 

First derivative El’R line shapes are efficiently calculated by numerical 

33 34 

diagonal izatlon nut hod. * 
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Figure ^ 
Figure 1 


Figure 2 


Figure 3 


Figure 4 


Ion h 

Graph of rotational correlation time t versus ordering parameter 
S for the breakdown of motlonally narrowed formulae. (— ) is 

breakdown for line separations; ( ) ts breakdown for line widths. 

Calculations wore made for lsotioplc Brownian rotation "ith simple 
one-term Maler-Saupe restoring potential. VOAA magnetic parameters 
were used. 

A comparison of simulated and experimental spectra for VOAA in 

m: PC at 79*C. ( ) is experimental spectrum. ( — ' Is calculated 

for isotropic Brownian rotation witli t * 4.1 x 10 sec/rad and 

S ■ -0.13. (••••) is calculated for Brownian rotation and aniso- 

9 

tropic viscosity with t • 4.1 x 10 sec/rad, Tj| *5.1 x 10 ^ 

sec /rad, and S * -0.13. 
z 

Comparison of simulated and experimental spectra for VOAA in B!'PC 

at 79°C and 38®C. ( ) are experimental spectra. ( ) are 

calculated with c' • 1.5, c' ■ 4. (••••) are calculated with 

sec psec 

c' - 1.5, c' - 10. At 79 P C, t - 4.1 x 10" 10 sec/rad, S » -0.13 

sec psec z 

at 38“C, t ■ 1.7 x 10 9 sec/rad, ■ -0.21. 

Comparison of simulated and experimental spectra for VCAA In Phase 

V at 69®C and 0®C. ( ) are experimental spectra. (-- ) are 

calculated with c’ «■ 1.5, e' "4. (••••) are calculated wit!; 

sec psec 

£' - 1.5, c' - 10. At 69°C, T - 2.6 x 10 _10 s/rad,3 - -0. U3; 

ovC pbuc 2 

at 0°C, T - 5.9 x 10" 9 sec/rad, S - -0.22. 

z 
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Figure 5. Graph «> 1 ordering parameter versus temperature T tor VOAA In 

V (• ’ 10) 0 - C' • 4 UPC (A - i ' - 10) 

puec psec pace 

A - ' *4). Sol Id lino la drawn lor Phase V and daubed lino 

peoc 

la for UK PC. » ' -1.5. 

sec 


Figure 6. Graph of rotational correlation time t (sec/rad) versus - lor 

VOAA in Phase V (s - t ' - 10; - ' - 4) and UKPC (A - i ' 

pace psec psec 

A - r' “4). Solid line la drawn for Phase V and dashed line is 

pHOC 

for BKPC. r ' -1.5. 
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